Band-filling effects on electron-phonon properties of normal and superconducting state 
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We address the effect of band filling on the effective electron mass m* and the superconducting 
critical temperature T c in a electron-phonon system. We compare the vertex corrected theory with 
the non-crossing approximation of the Holstein model within a local approximation. We identify 
two regions of the electron density where m* and T c are enhanced or decreased by the inclusion of 
the vertex diagrams. We show that the crossover between the enhancement at low density and the 
decrease towards half filling is almost independent of the microscopic electron-phonon parameters. 
These different behaviors are explained in terms of the net sign of the vertex diagrams which is 
positive at low densities and negative close to half filling. Predictions of the present theory for 
doped MgB2 , which is argued to be in the low density regime, are discussed. 



I. INTRODUCTION 



Recently, the discovery of superconductivity at T c ~ 39 K in magnesium diboride MgB2, 1 has stimulated a strong 
interest on the role of band-filling in phonon-driven superconductivity when the charge carriers stem from low electron 
or hole density in the conduction band. Classical examples of low density superconductors are doped semiconductors 
like GeTe which display critical temperatures T c of the order of 1 K or less. Such small T c values are understood as 
basically due to low values of the density of states at the Fermi level Nq which goes to zero as No cx y/E-p, where 
Ep is the Fermi energy. In this case, the electron-phonon coupling, which scales with the density of states, A = g 2 No 
where g is the electron-phonon matrix element, also vanishes as A cx \/Ep. This situation changes drastically for 
systems having a marked two-dimensional character of the conduction band at the Fermi level. In fact, an ideal 
two-dimensional system has a finite value of Nq even when Ep — > 0, leading in principle to a sizeable A also for very 
low charge carrier concentrations. 

Magnesium diboride belongs in principle to such class of materials. In fact, the peculiarity of MgB2 resides on its 
electronic band structure which is built up by three-dimensional 7r bands and quasi two-dimensional a bands. 2 ' 3 Both 
sets of bands have large bandwidth of the order of 5-10 eV, however the Fermi energy cuts the a band quite near 
the band edge leading to E-p of only 0.4-0.6 eV. Since the Fermi level crosses the a bands around the T point, the 
low energy physics of the a bands is essentially captured by a simple parabolic band structure corrected by a small 
interplane hopping element. 3 The electron-phonon coupling has been estimated by several groups and the common 
agreement is that the a bands are strongly coupled to lattice vibrations of E^ g symmetry and frequency o-> p h of order 
70 meV, while the it bands have a somewhat lower interaction. 4-6 Hence, the main features of MgB2 can be described 
by a quasi-two dimensional parabolic band with low density of holes interacting with the lattice vibrations. 

After the discovery of T c = 39 K in MgB2, the idea of varying the hole density of the a bands to possibly enhance 
T c has been the object of several works, both theoretical and experimental. A possible route is to provide additional 
charge, electrons or holes, to the system by chemical substitution of Mg or B atoms, although the rigid band scheme 
is somehow questionable. For instance, by substituting Mg with Al it is possible to further decrease the Fermi energy 
of the pure MgB2 7 ' 8 and consequently to explore the non adiabatic regime. It is interesting to see that by increasing 
the Al content the critical temperature decresase 7,8 but not at the rate expected by first principles calculations carried 
on with Eliashberg scheme 9 . 

An alternative way to chemical substitution is to search for related compounds with similar electronic properties. For 
example, insulating LiBC compounds have been theoretically predicted to display high values of the superconducting 
critical temperature T c for low hole-dopings corresponding to Ep < 0.5 eV. 10 

Such small values of Ep together with sizeable electron-phonon couplings define a highly non trivial situation 
for which the classical Migdal-Eliashberg (ME) theory becomes questionable. A major inadequacy of the classical 
theory is the breakdown of the adiabatic hypothesis occurring when the Fermi energy becomes of the same order or 
smaller than the characteristic phonon frequency scale Ep < Wph, 11 or when A becomes strong enough to give rise to 
polaronic effects. 12 So far this issue has been often associated with narrow band systems, as fullerenes, 13 cuprates 14 
and manganites, 15 where the Fermi energy is constrained by the small electronic bandwidth. 
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As discussed above, however, the nonadiabatic problem can be relevant also in large band materials with low carrier 
density provided a non vanishing electron-phonon coupling A. The nonadiabatic degree Lu p h/E-p can be thus tuned by 
varying the charge carrier density or alternatively the chemical potential through the electronic band up to the top 
(or the bottom) of the electronic states. 

Aim of the present paper is to investigate in details the onset of nonadiabatic effects in large band systems as 
function of the charge density and their implications on the superconducting and normal state properties. We find 
two regions of the electronic filling in which nonadiabatic contributions behave in qualitatively different manners: 
at sufficiently low density they enhance the electron-phonon effects in normal and superconducting state properties, 
while close to half-filling nonadiabatic effects basically lead to a weakening. 



II. NONADIABATIC THEORY: NON-CROSSING AND VERTEX CORRECTION APPROXIMATIONS 

In this section we present a brief summary of the theory of nonadiabatic superconductivity in term of the Green's 
function formalism 16 where the nonadiabatic effects will be introduced at different stages of complexity. An explicit 
analytical derivation of the superconducting equations will be provided within the local approximation. 

Let us consider a single band electron system interacting with phonons. In the normal state, the electron propagator 

is: 

G(k)- 1 = iw n -e k + e -S(fc), (1) 

where cu n = (2n + l)irT is the fcrmionic Matsubara frequency, T is the temperature, ek is the electron dispersion, 
eo the chemical potential and S(fe) is the electronic self-energy. In Eq. (I), the index k stands for the 4-momcntum 
(k, uj n ). In the following, we assume that the Coulomb contribution to the normal state self-energy has already been 
absorbed in the electron dispersion and the chemical potential. However, for the moment, we do not need to specify 
the form of e^. The electron self-energy is determined by the electron-phonon interaction: 

E(k) = J2v^(k,k')G(k'), (2) 

k' 

where J2 k , = —(T/N)J2 m J2is.'- Up to the first order of the irreducible electron-phonon interaction kernel of the 
self-energy V^(k,k') is given by: 17 ' 18 

V N {k, k') = V(k, k') + V(k, k') ^ V(k, k")G(k")G(k" -k + k'). (3) 

k" 

In the above expression, the lowest order electron-phonon interaction V(k, k') is given by V(k, k') = g(k, \t') 2 D(k — k') 
where g(k, k') is the electron-phonon matrix element and D{k — k') is the phonon propagator. In the following, we 
assume a phenomenological model for D(k — k'), in the same spirit therefore of the Migdal-Eliashberg theory in 
which the phononic degrees of freedom are taken from experiments while the electron Green's function is calculated 
self-consistently. Hence, we shall not consider the effects of the phonon self-energy. These are however negligible at 
sufficiently low electron densities, that is the region of main interest to us. 

If the second term in Eq. (3) is neglected, the so-called non-crossing approximation (NCA) is recovered, which, in 
the limit of infinite bandwidth, is equivalent to the Migdal-Eliashberg approximation. 11,19 When the full expression 
of Eq. (3) is considered, the system is described within the vertex-corrected approximation (VCA), equivalent to the 
theory of nonadiabatic electron-phonon interaction discussed in Refs. 16-18. For later convenience, and to maintain 
a closer connection with the notation used in Ref. 16-18, we define the vertex function 

P(k, k') = v ( k , k")G{k")G{k" -k + k'), (4) 

k" 

so that Eq. (3) is simply 

V N (k,k') = V(k,k')[l + P{k,k% (5) 

The extensions of NCA and VCA to describe the superconducting transition are readily obtained by either in- 
troducing the Nambu formalism or by computing the pairing susceptibility. We adopt here the latter formalism so 
that the temperature T c at which the pairing susceptibility diverges is given by the solution of the following secular 
equation: 
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<KQ = J2v A (k,k')G(k')G(-k')cj>(k'), 
k> 



(6) 



where <j){k) is the off-diagonal self-energy and the kernel V\ is: 

V A (k, k') = V(k, k') + V(k, k 1 ) ^ V(k, k")G(k")G(k" -k + k') 

k" 

+V(k, k') v ( k ", k)G{-k")G{-k" + k-k') + J2 v ( k > k")V(k", k')G(k")G(k" -k- k 1 ). (7) 

k" k" 

The NCA equations are obtained by neglecting all the nonadiabatic corrections so that VA(k,k') is simply V(k,k'). 
In this limiting case the electron-phonon interaction enters in the same way in the interaction kernels of both the 
self-energy and the superconducting pairing: V A (k, k') = Vys(k 7 k') = V(k 7 k'). This is however not true in the VCA 
framework where in principle V A (k, k') ^ V^,{k^k'). In particular the nonadiabatic terms add considerable structure 
to the pairing kernel but symmetry considerations permit to simplify the structure of V A {k, k'). In fact, if the system 
has inversion symmetry V(k, k') = V(~k,—k') and G(— k) — G{k)* where G(k)* is the complex conjugate of G(k). 
By using these relations it is easy to see that the third term on the right hand side of Eq. (7) is the complex conjugate 
of the second term. Their sum is therefore twice the real part of, for example, the second term. Consider now the 
last term of Eq. (7). If we replace the summed index k" with ~k" + k — k 1 then the resulting expression is equal to 
the complex conjugate of the original term. Then also the last term of the gap equation is real. These considerations 
permit to rewrite Eq. (7) in the more condensed form: 

V A (k, k') = V(k, k')[l + 2P R (k, k')] + C(k, k'), (8) 

where Pr(/c, k') is the real part of the vertex function P(fc, k') defined in Eq.(4), and 

C(k, k 1 ) = ^ V{k, k")V(k", k')G{k")G{k" -k- k'), (9) 

k" 

is the so-called cross correction, which is a real function of k and k' . Finally, reminding that G(k')G(—k') = |G(fc')| 2 , 
we can see that the kernel of the linearized gap equation in Eq. (6) is also real and without loss of generality we can 
choose 4>(k) as a real function of k. 

In this article, we study both NCA and VCA as a function of band filling for a model density of states (DOS) which 
is constant through the entire electron bandwidth D. When the chemical potential eo is close to the top (bottom) 
of the band, this model effectively mimics a hole (electron) doped two-dimensional parabolic band, provided w p h is 
smaller than the entire bandwidth. This modelization is thus a good starting point to eventually discuss filling effects 
in MgB2. Within the NCA framework, Marsiglio has already considered this model in his study of the superconducting 
instability. 20 What he found was that the critical temperature T c decreases as the chemical potential eo approaches 
the bottom or the top of the band. Naively, this reduction is expected since, close to the band edges, the average of 
the DOS over an energy window of order cj p h gets reduced. Quantitatively, however, the reduction of T c is affected 
also by non-rigid band effects induced by the electron-phonon coupling. 

Within the theory of nonadiabatic superconductivity, this situation is drastically modified by the non-trivial band- 
filling dependence of the nonadiabatic corrections. To illustrate this point, let us consider the vertex function P(k, k') 
defined in Eq. (4). In previous studies we have shown that, at half filling, the vertex correction has a marked 
dependence on the momentum transfer q = k k', being positive for i>F|q| < w p h and negative for «F|q| > w p h, where 
vf is the Fermi velocity. 16-18 If there is no preference to forward scatterings, the average over |q| of the vertex function 
is generally small and negative, reducing therefore the electron-phonon pairing and consequently T c as compared to 
NCA. 21 ' 22 As the chemical potential is moved towards the top or bottom of the band, the q-dependence is more 
and more weakened and, when the distance between the band edge and eo is sufficiently small, the vertex function 
assumes an overall positive sign with negligible momentum dependence. 23 In this regime therefore the nonadiabatic 
corrections tend to enhance the pairing, and the reduction of T c given by NCA is counter balanced by the nonadiabatic 
corrections. Close to the top or bottom of the band, we expect therefore a competition between band edge effects 
which depress T c and nonadiabatic corrections which enhance T c . Which of these two mechanisms prevails over the 
other depends on the chemical potential, on the bandwidth and on the electron-phonon interaction strength. It is 
matter of this article to provide quantitative results of this interesting problem. 

A. Local approximation 

To focus on the role of band-filling in enhancing or decreasing T c , we employ some approximations both in the 
description of the electron-phonon interaction and in the way of dealing with the normal and off-diagonal self-energies. 
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First, we assume that phonon dynamics is well described by a single Einstein frequency u> , so that the phonon 
propagator is 

D(q) = D(l) = --J%, (10) 

where lui — 2lirT. In addition, we consider the electron-phonon matrix element as independent of the momenta, 
i. e., g(k, k') = g. This latter assumption disregards situations in which the electron-phonon interaction is strongly 
momentum dependent such as in strongly correlated metals in which large momentum transfers are suppressed by 
the rigidity of the correlated electronic wave- function. 14 > 24 > 25 Previous studies have shown that in this situation the 
surviving small momentum scattering has deep influences on the nonadiabatic corrections. 16 In order to focus on 
the role of the band-filling effects, however, we assume in the following the electron-phonon matrix elements and the 
electronic Green's function to be independent of the electronic momenta. This assumption, which corresponds to 
the local approximation whose validity and limitations will be discussed later, permits the comparison with an exact 
solution in the limit of zero density. We note however that electronic models interacting with purely local interactions 
are widely studied as basis for the recent developed dynamical mean-field theory which provides in some cases an 
exact solution in infinite dimensions. 26-28 The local approximation can be viewed as an average out of the momentum 
dependence over the whole Brillouin zone of electron-phonon matrix elements and of the electronic propagators. Thus 
we formally replace the electron propagators G(k) by their local expressions G\ oc (n) = ^2^G(k,n) which depend 
only upon the Matsubara frequencies u n . The electron self-energy becomes also momentum independent and can be 
rewritten as: 



£(fc)->E(n) = iu n -iW(n)+ X (n), (11) 
where W(— n) = — W(n) and \{~ n ) = x{ n )- With this notation, G\ oc reduces to: 

GUn) = [ de N ® — , (12) 

J iW(n) -e + e Q - XW 

where N(e) is the bare electronic DOS and the integral over the energies e is performed over the whole bandwidth D. 
We consider here a constant bare DOS system N(e) = 1/D [-D/2 < e < D/2}. The value of the chemical potential 
£o = corresponds thus to a half-filling case whereas, for zero electron-phonon couplings, e = —D/2 and e = D/2 
represent respectively a completely empty and a completely filled band. When |cq — D/2\ <C D, this approximation 
simulates a quasi-two dimensional parabolic band. By employing these assumptions, Eq. (12) reduces to 

Gioc(n) = ^ f D/2 de-— L — = -hif(n) + g(n)}, (13) 

DJ-d/2 iW (n) - e + e - x{n) D 



where we have defined 



9{n) = \ In 



W(n)i + lD/2-e + x(n)} 2 } 



W(n) 2 + [D/2 + e - X (n)} 2 



By substituting Eq. (13) into the definition of Vn of Eq. (3), the resulting interaction depends now only on 
Matsubara frequencies and it reduces to: 

V N (k,k') -» V N (n,m) = g 2 D(n - m)[l + AP(n, m)], (16) 

where A = g 2 N(0) = g 2 / D is the electron-phonon coupling constant and P(n,m) is the vertex function: 

P{n, m) = -DTj2 D ( n - l)G loc (l)G loc (l -n + m) = P R (n, m) + iPi{n, m), (17) 



where 



P R (n, m ) = -^J2 D ( n - l)[9(l)g(l - n + m) - f(l)f(l -n + m)} (18) 
l 
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and 



Pi(n, m) = - J Yl D ( n ~ l )\f(W - n + m) + f(l - n + m)g(l)} (19) 
l 

are the real and imaginary parts of the vertex function. Note that, since g{— n) = g(n) and /(— n) = —f(n), Eqs. 
(18), (19) imply Pr(— n, — to) = Pr(ti, m) and Pi(— n, — m) = —Pi(n,m). Now, by applying the local approximation 
to Eq. (2) and Eq. (11) we can write in a closed form the integral equations for W(n) and x( n ) : 

W{n) = u) n - \T^D(n - m){[l + AP R (n, m)]/(m) + \Pi(n, m)g(m)}, (20) 

m 

X (n) = \T^D{n - m){[l + AP R (n, m)]g{m) - APi(n, m)/(m)}. (21) 

m 

The above equations are completed by the relation connecting the electron occupation number n with the chemical 
potential e : 

» = 1-0&W (22) 

n 

where the factor 2 stems from the spin degeneracy. It is important to notice that, form Eq.(22), half-filling (n = 1) 
is achieved by g(n) — 0, which implies that the imaginary part of the vertex function Pi, Eq. (19), is zero. This is 
consistent with x( n ) — an d e o = [see Eqs. (15), (21)]. 

Finally we consider the local approximation in the superconductive instability equations, Eqs. (6), (7). By follow- 
ing the previous steps, and by using Eqs.(13)-(15), it is straightforward to deduce that the gap equation becomes 
momentum independent and that it reduces to: 

4>{n) = - T J2v A (n,m)f(m)^- (23) 

where the pairing interaction in the local approximation reads 

V A {n,m) = X[l + 2XP n (n,m)]D(n - m) + X 2 C{n,m). (24) 

In the above expression, P R is the real part of the vertex function given in Eq. (18), while C(n, m) is the nonadiabatic 
cross correction resulting from the local approximation of the last term of Eq. (7): 



C(n, m) = -TDj2 D ( n - l ) D ( l - m)G loc (l)Gi oc (l — n — to 



= ~ Y, D ( n - W - m)[g{l)g{l - n - m) - f(l)f(l -n- to)]. (25) 
i 



When Eqs. (14, 15, 18, 19, 25) are plugged into Eqs (20-23), we obtain a closed set of equations which can be numerically 
solved for both NCA and VCA to obtain the normal state self-energy and the superconducting transition. The 
numerical solution is achieved by using usual procedures: first the normal state quantities W(n), xi n ) are obtained 
by iteration by constraining the electron density n to have a given value; in a second step the gap equation, Eq. (23), 
is rewritten an an eigenvalue equation and T c is obtained when the highest eigenvalue becomes equal to the unity. 

Let us stress that this theory could provide a route to analyze the data from experiments via a theory which is 
valid even in a very low filling regime where Eliashberg and Mc Millan formula for T c are no longer valid. In fact 
given an clcctron-phonon coupling g(k, k'), as well as a calculated phonon spectrum, it is possible to evaluate an 
effective local interaction Va by averaging over the entire Brillouin zone the product V(k,k') = <?(k, k') 2 £>(fc — k') 
where D(k — k') is the corresponding phonon propagator. Notice that in this way filling effects will enter both in 
filling dependence of el-ph coupling and of the renormalized phonon frequencies. In this paper we conversely focus on 
the explicit band-filling dependence of the Fermi energy which are the driving effects at low density. To this aim wc 
therefore choose a simple lorentzian form for V(k, k'). 
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FIG. 1. Electron density as a function of the chemical potential eo evaluated according NCA and the VCA approximation 
for finite (ooo/D = 0.1) and almost zero (luq/D = 0.01) adiabatic parameter. Empty circles represent the exact result of the 
zero density limit evaluated by means of the continued fraction within the local theory. Electron-phonon coupling is set A = 1. 



III. NORMAL STATE PROPERTIES 



Our main aim is to investigate how the electron-phonon properties (effective mass, superconducting pairing, etc. . . ) 
depend, in both NCA and VCA, on varying the chemical potential as function of the electronic filling in comparison to 
the phonon energy scale luq and to the large electron energy scale provided by the total bandwidth D. In this section 
we focus on normal state properties while the superconducting pairing will be discussed in Section IV. In relation to 
both these quantities we shall distinguish two regimes: one at low electron filling n where vertex diagrams enhance 
the effective electron-phonon coupling associated to the renormalization of the electron mass and the particle-particle 
superconducting pairing; one at large n where the renormalization due to the vertex processes leads to an effective 
reduction of the coupling. We shall see that, while the magnitude of these electron-phonon effects depends on the on 
the strength of A and on the size of the adiabatic parameter ojo/D, the typical filling where the crossover between 
these two regimes occurs is only weakly dependent on the electron-phonon properties, suggesting that it is mainly 
determined by purely electronic reasons. The relation of these results with the sign and magnitude of the vertex 
function will be discussed. 

Concerning the normal state we present results for two physical quantities: the electron density as a function of the 
chemical potential and the effective mass ratio. The first quantity is obtained through Eq. (22) while the effective 
mass ratio is evaluated as low energy limit of the rcnormalized electronic frequencies [Eq. (20)] : 

m* W(n) 
m uj n 

This expression provides a good estimation to the band mass as far as the temperature is sufficiently low. We have 
iteratively solved equations (20), (21) at several temperatures until we have approached the asymptotic low temperature 
mass ratio. In this case we found that also the electron density is converged to its low temperature limit. Calculations 
have been carried out at several values of the chemical potential for very low {wq/D — 0.01) and moderately large 
(loo/D = 0.1) values of the phonon frequency and for two values of the coupling constant representative of weak 
(A = 0.5) and strong (A = 1.0) coupling. 

In figure 1 we report the behaviour of the electron density as a function of the chemical potential evaluated using 
both NCA and VCA. No significant difference between the two approximations is observed. Note that, for small 
values of n, e falls below the lower edge e = —D/2 of the bare conduction band. This is especially evident when 
u>o/D = 0.1 for which eo < —D/2 already for n < 0.15. In the zero density limit the chemical potential approaches 
the bottom band edge which corresponds to the ground state energy for a single electron in interaction with phonons. 
This energy is lowered with respect to the bare value (eo/D = —0.5) by the electron-phonon interaction. Note that 
in the n = limit of this local theory the electron and phonon properties can be evaluated in a controlled way as 
well as in NCA and VCA by using the continued fraction technique (NCA and VCA correspond to different specific 
approximations of the continued fraction). 27 This is quite similar to the d = oo case where the local ansatz is enforced 
by the infinite dimension limit. 26 The zero density exact result of the continued fraction is also reported in the figure 
(circles). Note that, at least for these value of ujo/D and A the vertex corrected theory provides already a rather good 
estimate of the exact results at zero density. 29 We can expect thus a polaron crossover around some critical value 
of the electron-phonon coupling A c which depends on the DOS shape. This feature however can not be reproduced 
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FIG. 2. Effective mass ratio for u>o/D — 0.01 according the non-crossing (NCA) and the vertex corrected (VCA) approxima- 
tions. In the zero electron density limit data are also shown for exact calculations by means of continued fraction and for the 
second order perturbation theory (2 pt). 



in NCA and VCA since it implies resummation of higher (infinite) order vertex diagrams. 27 The reliability of our 
approach is thus limited to the metallic regime of the electron-phonon system here considered. The value A c ~ 1.544, 
corresponding to the polaron transition in the adiabatic one particle electron-phonon system with rectangular DOS, 
represents thus an evaluation of upper limit of validity of our analysis near the adiabatic limit. Note however that 
for A < A c there is no reduction of the zero density ground state energy due to the electron-phonon interaction in the 
strictly adiabatic limit i.e. when loq/D — > 30 . This result, which includes all order vertex corrections, holds true even 
when the local approximation is relaxed. Therefore we expect that frozen phonon approximation cannot reproduce 
this low density behaviour at least for couplings outside the polaronic regime. 

After having discussed the n-eo phase diagram and the range of validity of our analysis, we address now the filling 
dependence of the effective mass ratio m*/m. In figures 2 and 3 the results of the calculations for ojq/D = 0.01 
and loq/D = 0.1 are respectively reported. Because of band edge effects the effective electron mass m* is reduced 
in NCA as the filling approaches n — > 0. Increasing the phonon frequency leads to a broadening of this trend. The 
introduction of vertex diagrams can affect this scenario. In particular, as shown in Fig. 3 (lower panel), for large 
enough electron-phonon coupling and adiabatic parameter loq / D the effective electron mass m* in VCA can be larger 
at the band bottom than at half-filling suggesting an increasing of the effective electron-phonon coupling in this regime 
as due to the vertex diagrams. 

The comparison between NCA and VCA can provide an estimate of the vertex effects alone. In particular for both 
low and intermediate phonon frequencies we can distinguish a low doping regime in which the inclusion of the vertex 
diagrams increases the effective mass 29 and a high doping regime in which effective mass is reduced by the inclusion 
of the vertex corrections. 33 

In order to single out the effects of the nonadiabatic diagrams as function of the electron filling, we plot in Fig. 
4(a) the relative differences for VCA and NCA results concerning the renormalization of the electronic mass m* 
[(m* VCA — m* NCA )/m* NCA ]. This plot shows clearly the two regimes previously mentioned: a low electron density 
regime where the introduction of nonadiabatic vertex diagrams (VCA) leads to an enhancement of the electron-phonon 
coupling with respect to NCA for what concerns single particle properties (m*/m); and a large electron density regime 
where the opposite holds true. An interesting results is that the crossover filling n c ~ 0.15 — 0.2 where nonadiabatic 
vertex effects change from positive [(m* VCA - m* NCA )/m* NCA > 0] to negative [(m* VCA - m* NCA )/m* NCA < 0] is 
almost independent of the bare electron-phonon coupling A and of the adiabatic ratio ujq/D. We would like to stress 
again that for these values of n the chemical potential eo, in particular for ujq = 0.01D, can be quite far from the 
bottom of the band with respect to the phonon energy, so that the crossover is not simply related to band edge effects 
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FIG. 3. Effective mass ratio for ojq/D = 0.1. Notations are the same of Fig. 2. 
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the electron filling. Panel (b): sign of the real part of vertex function Pn(n,m) [Eq. (18)] evaluated at zero frequencies n = 0, 
m = 0. 
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as described by a rigid band picture. 

In order to better clarify the origin of this crossover from a VCA m*-enhancement region (n < 0.15 — 0.2) to a VCA 
m*-reduction one (n > 0.15 — 0.2), we plot in Fig. 4(b) the size of the vertex function Pr(ti, m) [Eq. (18)] evaluated 
in the static limit n = 0, m = 0. Previous studies showed that the static limit can be indeed representative of the 
effects of the vertex function on static quantities, like m*. The comparison of panels (a) and (b) of Fig. 4 points out 
striking similarities between the filling behaviour of the two quantities. It is thus easy to identify the sign of Pr as 
the direct origin of the enhancement or suppression of the VCA electronic mass m* with respect to NCA. Once more 
we would like to stress that while the magnitude of the discrepancy between NCA and VCA (panel a) or Pr (panel b) 
depend on the adiabatic parameter ujq/D and on the electron-phonon strength A, the electron filling n c ~ 0.15 — 0.2 
where vertex diagrams are almost perfectly balanced is only weakly dependent on electron-phonon properties. 

An other interesting issue to be investigated is the dependence of the electron-phonon renormalization on the degree 
of the diagrammatic approximations employed. Such a dependence can be traced out in a controlled way in the zero 
electron density limit where different theories correspond to different degrees of diagrammatic approximations. 27,29 
We compare thus the NCA and VCA with: i) the second order perturbation theory (2 pt) where only the lowest order 
diagram with bare propagators is taken into account, namely in Eqs. (20)-(21) the r.h.s. is evaluated with W n = oj n , 
X = and Pr = Pj = 0; ii) and with the exact result for the local theory obtained by the infinite resummation of the 
continued fraction. 

For all the adiabatic parameters ujo/D and coupling A here considered (Figs. 2, 3) we observe a monotonous increase 
of the effective mass enhancement m*/m from the second order perturbation theory, to NCA, VCA, ending up with the 
exact result. This suggests that vertex diagrams in the zero density limit are positive defined (the well defined positive 
sign of the vertex function in the zero density limit was analytically pointed out in Rcf. 34 ) and the resummation 
of higher and higher vertex diagrams, as taken into account by the continued fraction, converges uniformly for the 
values of A and loq/D here considered. In this situation the vertex corrected theory can be considered qualitatively 
representative of the infinite resummation exact result. 29 This picture is expected to be less and less accurate as soon 
as one approaches the polaronic regime. For instance we note that in the almost adiabatic case (ojo/D = 0.01) VCA 
provides already a good estimation of the zero density exact result even for large coupling A = 1.0. This can be 
understood by reminding that electron-phonon polaron correlations are frozen in the adiabatic limit up to A < A c , 
while they diverge at A = A c 27,30 On the contrary for loq/D = 0.1 the approaching of the polaron crossover is reflected 
in different properties (here for instance in the electron mass renormalization) even for A < A c . This means that the 
qualitative trend from 2 pt to NCA, VCA and to the exact result at n = is still preserved, but it lacks under the 
quantitative side as far as A get closer to A c . 29 

Summarizing the results of this section, we have identified two distinct regions of the electronic filling where the 
vertex diagrams are respectively positive and negative, leading to a corresponding increasing or reduction of m* with 
respect to NCA. We have also shown that the extension of this filling region does not vanish in the adiabatic limit. The 
positivity of the vertex diagrams can be better understood in the zero density limit. This is a quite general result and 
we believe that the physics of this limit can be qualitatively representative of the low doping region (n < 0.15 — 0.2) at 
least in the case where no phonon renormalization effects are not explicitly taken into account (note however that the 
phonon screening vanishes as the charge density goes to zero). Once again we stress that at large filling n > 0.15 — 0.2 
the system behaves in a quite different way, since for instance the inclusion of vertex diagrams leads to a reduction of 
the effective mass m* instead of an enhancement of it. 



IV. CRITICAL TEMPERATURE 

The competition between band edge effects and nonadiabatic contributions on the critical temperature shares quite 
similarities with the case of the effective electron mass m*. In Figures 5 and 6 we report the results for T c in both 
NCA and VCA frameworks, respectively for w /D = 0.01 and for ujq/D = 0.1, as function of the electron density n. 

In analogy with the renormalized electron mass m*/m, the NCA results show that T c reduces as the electron 
occupation number n deviates from half-filling (n = 1). This can be again naively understood in terms of reduced 
density of states as the band edge is approached within a phonon energy window. The resemblance with the n- 
dependence of m* /m is maintained also when nonadiabatic corrections are included: at half filling T c is reduced with 
respect to NCA while, at sufficiently low n, T c is enhanced. This result is in agreement with previous studies on the 
density effects on nonadiabatic corrections. 21,23 At a first sight it could appear surprising that an enhancement of the 
critical temperature is found in the presence of a corresponding enhancement of the effective mass m* since strong 
renormalization of m* is usually associated with a depletion of spectral weight. However it was shown that even in 
the adiabatic theory in the strong coupling limit both critical temperature (T c cx V~X) and effective mass [m* jm cx A) 
increase. 31 In addition it was pointed out in Ref. 32 and later explicitly shown in Refs. 16-18 , that a different mechanism 
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FIG. 5. Critical temperature T c /uo vs. the electron filling n for loq/D — 0.01 in NCA and VCA. 
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based on nonadiabatic contributions leads to a reduction of the spectral weight which is counterbalanced by the 
additional opening of nonadiabatic interaction channels in the Cooper paring even at intermediate coupling. 

Although m* /m and T c behave in a qualitative similar manner, interesting differences can be singled out by looking 
at the relative change of T c due to the vertex diagram inclusion: (T^ CA — T^ CA ) /T^ CA (Fig. 7a). We observe that 
the filling n c at which the inclusion of vertex diagrams changes its character from positive to negative with respect to 
the superconducting critical temperature is now quite sensitive of the clcctron-phonon coupling A and of the adiabatic 
ratio u>o/D. In particular n c increases as A and ujq/D increase. This is in contrast with the case of the electron mass 
where n c was almost independent of electron-phonon quantities. There are two basilar reasons which account for this 
difference. First, we remind that, while the electron mass renormalization is in good approximation a static quantity 
(e.g. it is evaluated at zero frequency), the calculation of T c involves the knowledge of the electron-phonon kernel 
(and of the vertex function) in an energy window ~ ±wo- The overall structure of the vertex function can be thus 
remarkably dependent on A and ujq/D affecting the global balance for which T^ GA — T C NCA . An additional source of 
discrepancy comes from the fact that T c is by definition a finite temperature quantity. The value of T c itself is thus 
a X-uj / D dependent quantities and finite temperature effects can be drastically important in Fig. 7a. 

In order to single out these finite temperature effects we evaluate the maximum eigenvalue w max of the supercon- 
ducting electron-phonon kernel Vj\(n, m) [Eq. (24)] at low temperature in the normal state. This provides thus an 
estimate of the superconducting pairing which is not affected by finite temperature effects. A zoom of the interesting 
region [(v^ A - v^)/v%££ ~ 0] is shown in Fig. 7c as well as a corresponding zoom of (T C VCA - T C NCA )/T C NCA - 0. 
As a first remark we note that, contrary to Fig. 7b, in Fig. 7c there is only a weak dependence of n c on the electron- 
phonon coupling A. As matter of fact the dependence on A in Fig. 7b mainly stems from finite temperature effects, 
with a large A = 1.0 associated with higher temperatures (T c ~ 0.1 — 0.12) in comparison of T c ~ 0.03 — 0.04 for 
small A = 0.5. Once disregarded finite temperature effects, we still note a significant dependence on the adiabatic 
ratio u>o/D in Fig. 7c which does not appear for instance in Fig. 4. As above discussed, we trace back this residual 
dependence on ojo/D to the important role of the frequency dependence in the superconducting kernel. Interesting, 
in the almost adiabatic regime ujq/D = 0.01, the same value of n c ~ 0.15 as found for the electronic mass and for 
Pr{h> = 0) m = 0) is recovered, signalizing the static nature of that value. It is significant that both taking into 
account the retarded nature of the superconducting pairing and finite temperature effects yield an increase of the low 
doping region where nonadiabatic vertex diagram enhance the electron-phonon coupling. The value n < 0.15 — 0.2 
can be thus considered as a conservative estimate of this regime. 

In order to focus on the effects of nonadiabatic electron-phonon contributions, we have not considered in this work 
the role of Coulomb repulsion between electrons. In the simplest scheme, the electron-electron repulsion can be 
included by adding an Hubbard U term in Eq.(6), as done in Ref. 20. By treating U as a parameter independent of 
electron filling we have performed calculations of T c in both NCA and VCA. We have found that the inclusion of U 
leads to an overall reduction of T c and to a magnification of the differences between NCA and VCA reported in this 
section. However, no qualitatively new physics is induced by U, although we expect an additional reduction of T c 
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when the electron density is sufficiently low to poorly screen the Coulomb interaction. However, the role of screening 
was not the aim of this paper, and its effects will be investigated in a future publication. 

V. DISCUSSION AND CONCLUSIONS 

In this paper we have addressed the role of band filling in some normal and superconducting properties of electron- 
phonon systems. In particular, within a local approximation, we have identified two regions of the electron band 
filling in which electron-phonon properties behave in qualitatively different manners. At sufficiently low density, the 
inclusion of vertex diagrams beyond the non-crossing approximation leads to a strengthening of the electron-phonon 
effects reflected by an enhancement of the effective electron mass m* and of the superconducting critical temperature 
T c . On the contrary, close to half filling, the vertex corrected approximation results in a net decrease of m* and T c 
with respect to the non-crossing approximation. We have shown that these different behaviors are governed by the 
net sign of the vertex diagrams which is positive in the low density regime and negative close to half-filling. The 
density value n c ~ 0.15 — 0.2 separating these two regimes depends only weakly on the microscopic electron-phonon 
parameters which, on the other hand, rule the size of the corresponding enhancement or decrease of m* and T c . 

The results for the low electron density regime are compared with the zero density limit for which the exact solution 
of this local problem is known. In this limit we have shown that the inclusion of diagrams with higher and higher 
order, from second order perturbation theory to the non-crossing approximation and to the vertex corrected theory, 
monotonically approaches the exact results. Such an improvement as function of the degree of approximation suggests 
that the positive sign of the vertex function leads to a series with positive terms, which should steadily approach the 
exact solution if a resummation of infinite order diagrams were performed. 29 

The enhancement of the effective electronic mass m* and of the superconducting critical temperature T c in the low 
electron density region is expected to be relevant in twodimensional systems where small charge density is accompanied 
by a finite density of states at the bottom (or top) of the band. Thus, in principle, our results could concern the 
high-T c compound MgB2 and the hypothetical doped LiBC recently proposed to superconduct with even higher T c . 10 
In particular in MgB 2 the hole carrier density related to the a bands n h ~ 0.08 < n c 3 sets the system in a region where 
these effects are relevant, although the large bandwidth of the a bands in MgB 2 compared with the relevant phonon 
frequency leads to ujq/D <~ 0.01, corresponding to a quite weak effect on m* and T c . These effects can be however 
further enhanced via Al substitution. Interestingly, the experimentally observed reduction of T c as function of Al 
doping 7,8 cannot be accounted within a McMillan scheme 9 even when the Coulomb pseudopotcntial [i* is introduced 
to account for T c ~ 39 K in the undoped MgB 2 compound. 

A quantitative enhancement of T c from vertex corrections in MgB 2 could be argued by considering two distinct 
effects not considered in the present paper namely: 

a) the intrinsic momentum dependence of the bare electron-phonon matrix elements; b) the momentum dependence 
of the renormalized electron-phonon vertex. 

Let us discuss them separately. 

a) Several first principle calculations have shown that the electron-phonon coupling in MgB 2 has a strong momentum 
dependence and is particularly strong near the T point at q = 0. 4 ' 5 In this situation, neglecting the momentum 
dependence of the electron-phonon matrix elements t/k.k' ~ 9 becomes a poor ansatz just as in the case of strongly 
correlated electron systems. ' ' A proper inclusion of the predominance of forward scattering at small q has shown 
to lead to an intrinsic selection of the positive sign for the vertex diagrams and to favor a corresponding enhancement 
ofm*andT c . 16 

b) Even neglecting the momentum dependence of the bare electron-phonon matrix elements the renormalized vertex 
acquires a momentum dependency due to vertex corrections. The local approximation treat these effects by averaging 
the intrinsic dependence of the vertex diagrams on the momentum transfer q = k k'. In extreme cases the local 
assumption can ill-estimate the total phase space available for these vertex processes and hence the magnitude of 
the vertex diagrams. This is the most evident for instance in the low density adiabatic regime where the available 
transferred momenta are limited, in weak coupling regime, by ~ |q| < 2fcp, while the local approximation would 
integrated them over the whole Brillouin zone. In this limit for instance it is easy to show that the local approximation 
provide an underestimate for the magnitude of the vertex diagrams. 

Both effects could be separately considered to account properly of vertex corrections in MgB 2 . Work in this direction 
is in progress. 

This work was partially supported by the projects COFIN-2001 MIUR and PRA-UMBRA INFM. We thank G. 
Profeta for useful discussions. 
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